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We propose a simple effective model to describe FeAs superconductors. This model is based on 
the assumption of a local spin-density-wave (SDW) order, with its magnetization direction allowed 
to fluctuate. It is shown that the long-range order with momentum Q = (tt, tt) is generally unstable 
in competing with the kinetic energy of the charge carriers. A true weak SDW order is formed 
in the undoped case with an additional momentum shift Qj, = (tt, O) due to the peculiar Fermi 
surface nesting. In the doped case, the fluctuating long-range order driven by kinetic energy can 
naturally result in a d-wave superconducting condensation. Such low-energy physics is protected 
by the presence of the local SDW which sustains some kind of "Mott gaps" for the multiband 
d-electrons near the Fermi energy. 
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Introduction. The recent discovery of the iron-based 
superconductors [l|, 0, [^, 0, 01 has stimulated a lot of 
interest concerning underlying mechanism for supercon- 
ductivity in this new superconductor family. It has been 
established by various measurcmcnts[^ 0, 01 that there 
exists an SDW order in the undoped LaOFeAs compound 
below T ~ 150 K, which quickly disappears with the 
electrons doped into the system, where the supercon- 
ducting phase starts to set in with transition tempera- 
tures being raised beyond 50 K.[^, [§] While the general 
phase diagram[l|, d, 0, |3j Q reminds us some interest- 
ing similarities with the cuprate superconductors, the d- 
electrons on Fe seem much more itinerant with multi- 
orbitals crossing the Fermi level as indicated by the band 
structure calculations, [lol 11, [13, 13, 14, [ij [l^ com- 
pared to the isolated Gu^d 2_ 2"02p^a antibonding or- 
bitals in the cupratcs. Many theoretical pro posals are 
based on itinerant approaches [TtI [isl . [Til. l20j with the 
emphasis on the important role played by various mag- 
netic fluctuations, while some conjectures are also made 
from the side of large-spin Mott insulators. (2l| 



The LDA calculations [ij, ll5[ have found an energet- 
ically robust SDW state at the antiferromagnetic (AF) 
momentum Q = (tt, tt) with a large Fe moment ~ 
per site, but experimentally only a weak SDW ordering 
with a different magnetic momentum Qs = (tt, 0), which 
further doubles the unit cell of the former SDW state 
composed of two Fe per cell, has been identified^, @1 in 
the undoped case. The latter SDW (called stripe type 
below) was predicted by the first principle band struc- 
ture calculation!^ due to the nesting Fermi surfaces of 
the hole and electron pockets, which is much more "frag- 
ile" and easily destroyed as the doped electrons fill up 
the small hole pockets at small doping. 

In this paper, we will make a very simple proposal by 
assuming that the SDW state with an AF momentum 
Q remains strong locally in both the undoped and small 
doped regime. The corresponding profile of the electron 
density of states is schematically illustrated in Fig. 1. 
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FIG. 1: (color online) Schematic proflle of the density of states 
near the Fermi energy in the presence of the local SDW mean- 
flelds Ma and Mf,, which decide the "charge gaps" for the a- 
orbital and 6-orbital bands according to Hq [(O]. The lower 
half of each band (the so-called a-band, see text) is assumed 
below the Fermi energy Ef and is flUed, while the upper half 
(the /3-band) is only partially flUed. Here a-bands refer to the 
bands passing the Fermi level, which are more dxy, dxz, dyz 
like, and &-bands are more of d^2_y2 character according to 
the LDA resuh.Hfil 



Based on this minimal model, we can show that the long- 
range part of this SDW is actually generically unstable, 
by coupling to the charge carriers near the Fermi level. 
Such instability can result in a weak stripe-type SDW 
order in the undoped case where the nesting Fermi sur- 
faces are present, and naturally a superconducting con- 
densation when the former is destroyed at small doping. 
In this whole regime, however, the local order of such an 
SDW remains robust to "protect" the low-energy physics, 
which resembles the role of the "Mott gap" in the cuprate 
superconductors, [i^l 

Minimal Model. In our model Hamiltonian Hcs = 



band 



Hi, the first term is a tight-binding model 



band 



ij,a,b,a 



(1) 



which describes the electron effective hoppings between 
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FIG. 2: (color online) (a) The iron atoms form a square lat- 
tice with the dashed diamond denoting the unit cell of the 
SDW order state with magnetic momentum Q. (b) The mag- 
netic Brillouin zone (BZ) is illustrated by the dashed diamond 
shape. Two circles at the F and M points are hole and elec- 
tron pockets connected by a momentum Qs. The smaller 
square is a reduced BZ corresponding to an enlarged unit cell 
indicated by the largest square in (a). 



the d-orbitals of the Fe ions on square lattice (Fig. 2(a)) 
at the nearest neighboring (NN) and next nearest neigh- 
boring (NNN) sites, including intra- and inter-orbital 
hoppings with the superscripts a, b specifying the or- 
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bitals. There have been several proposals 
for i?band based on the LDA calculations in order to cap- 
ture the relevant bands near the Fermi energy. 

The second term Hi reflects the effective influence 
of the Coulomb interaction on the d-electrons, which 
includes the on-site and NN repulsions as well as the 
Hund's rule ferromagnetic coupling, and has the follow- 
ing "mean-field" look 



Ft 



(2) 



where M^a denotes an effective SDW mean-field felt by 
the electron spin Sta of the orbital a at site i. In an SDW 
ordered state, one would have Mj^ cx X]fc(5!^a) '^'h' i^ib) + 
Ua (Sia), where J'^'' is the Hund's rule coupling constant 
and Ua is the on-site repulsion. We shall make the fol- 
lowing ansatz 



(3) 



where a single unit vector (— l)*ni will describe the true 
polarization direction of Mia- Here the staggered factor 
(— 1)' is introduced such that = n will correspond to a 
true AF order, but in a general case rij will not be fixed 
around a particular direction as only the relative change 
AaXii = iii+Q, — n; (a = X, y) will enter the Hamiltonian 
as shown below. 

The local SDW field Ala will be assumed large accord- 
ing to the LDA calculation. (13. 15 1 with the magneti- 
zation at different orbitals tightly aligned together by 



the Hund's coupling. Under this assumption, the long- 
wavelength, low-energy fluctuations of may be treated 
as an independent degree of freedom, will be self- 
consistently determined by coupling to the electrons near 
the Fermi level in 7?eff- 

Effective theory. One may redefine as the new z-axis 
for the spin index of the electron spinor operator: 

Cia = UiCL^a (4) 

with = (aIa,T'^L,i) an SU(2) rotation 



U]-n.i ■aU,=&, 



Then Hi simply reduces to 



Hi = -J2Mai-iysi 

i,a 



(5) 



(6) 



(7) 



while Hhand becomes 

We may further rewrite 

Hctf = Ho + Hi (8) 

where Hq = Hi + -ffband[t^/C/j = 1] is simply an SDW 
mean-field Hamiltonian for the multibands, which can 
be diagonalized as 



k,a 

J2 ^£ ("L«ka - /3L/3ka) + const. (9) 



k.a 



by a canonical transformation Oka = u^ctka ~ 

V^a^Pl^a, Ok+Qa = v'^O-^&Ua + u'^Pka with U^. = 

[(l-er/^k)/2]'^', < = [(l+er/i?iJ)/2]'/'. Here 



^k^^ (^k±^k+Q)/2and 



{en' + iMa/2f 



(10) 



where denotes the bare spectrum determined by ([T]) 
(setting the chemical potential n — 0). Note that the 
band label a here can be different from the original orbital 
label in ^ because of the mixture of orbitals, and in 
obtaining (fTO|) . the same Ma is assumed for the mixed 
orbitals. Here k is defined in the magnetic Brillouin zone 
(BZ) with the magnetic momentum Q = (tt, tt), which 
coincides with the BZ of two irons per unit cell (Fig. 
2(b)). 

Now let us consider the term with UjUj ^ 1: 



^1 = - E ^f^a {ujU, - l) 



(11) 
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Wc shall focus on the case in which the a-bands are all 
filled up by the electrons, and the Fermi energy is lo- 
cated in some of /3-bands, which corresponds to both the 
undoped and electron-doped (or slightly hole-doped) sit- 
uations as illustrated by Fig. 1. After integrating out 
the a-bands and by assuming Aciii is small, Hi may be 
simplified [23! (at large Ma) to 



i/i ^ i ^ Vsk • (Dq)„ 



9k+q + 0-6'k+q) ^k+q,cr/3k,-a- 
(12) 



k,q,CT 



where Dn is a Fourier transformation of 



D, 



-{An, X n, 



(13) 



and 



'k+ql^k+qj 



is the step function restricting k + q 
within (outside) the magnetic BZ. Note that for sim- 
plicity we have omitted the band indices in (jl2p . where 
the spiral field Dq will couple to the electrons from all 
/?-bands near the Fermi level. It shows that the twist 
of Ui is kinetic energy driven. Such a spiral twist is 
balanced by the "superexchange" term in Hi, where 

Jeff ~ l/^Eka(l - "ka)(^k")V^^k prOvidcS thc Spin 

stiffness against the twist of thc SDW order and is mainly 
contributed by the filled a-bands and is reduced with the 
increasing filling in the /3-bands. 

Some technical remarks with regarding the deriva- 
tion of are in order. Strictly speaking, the dy- 
namic field rii should be introduced in a path-integral 
formalism, 2^ 24 1 where the temporal term U|^tU^ will 
also enter the Lagrangian. As shown previously for the 
one-band Hubbard model, 23, 2^ the terms like ujdtUi 
will play an important role in determining thc dynam- 
ics of rii at half-filling where an SDW long-range or- 
der is present. In this case, by expanding n.; around 



iq(nq X z) • CT, the propagator of Dq reduces 



to I?''(q, <) ~ qq(nq (t) • riq (0)) which is proportional 
to the spin-wave propagator with the coupling to doped 
particles vanishing as q — > 0. However, once a finite den- 
sity of thc holes or electrons are doped into the system, 
the dynamic spiral fluctuation of Dq which couples to 
the doped charge carriers in ()12p will become dominant 
over the other fluctuation terms like ujdtUi, leading to 
intrinsic instabilities [2^ of the system (see below). 

The propagator V{q,t) = -i (TtTr [Dq(t)D_q(0)]) 
can be obtained by integrating out the /?-band 



P(q,c^) ~ -l/[n^(q. 



Jos] 



(14) 



where 11^ (q,^) is the "bubble" diagram contribution 
of the /3-band electrons, with thc "bare" 11^(0, 0) ~ 
—VpNp ~ —t where vp and Np are the Fermi velocity 
and density of states based on Hq , which is expected to 



dominate over Joff on general grounds. [2 5| reflecting the 
fact that the hopping energy is always dominant over the 
superexchange energy locally, leading to the instability 
of the AF long-range order. 

"Stripe" instability in the undoped case. According to 
the band structure calculations, IJ, |15[ in thc undoped 
case, there are hole-pockets around the F point and elec- 
tron pockets around the M point around the Fermi level, 
which are connected by a particular nesting momentum 
Qs, and can be described by Hq with the proper choice 
of Then, to take advantage of the enhanced response 
function at Q^, in Hq, the interaction term Hi will natu- 
rally induce 



Vek- (DqJ ^0 



(15) 



and result in a mean-field ^/3j^_|_q o-Z^k.-o-y 7^ with a 

\/2 X \/2 folding of the BZ as shown in Fig. 2(b). Con- 
sequently a gap will be opened up at the Fermi energy 
to stabilize P(q, w) in p4|) . In contrast to thc unstable 
SDW with the AF momentum Q, thc true magnetic or- 
dering (|15|1 is realized with an additional spiral twist at 
momentum Qs, which, known as a stripe- type SDW, was 
first predicted in Ref. Q| and has been recently confirmed 
experimentally by neutron scattering. 0, 0| 

Superconducting instability. With introducing a small 
amount of doped electrons, thc Fermi energy will move 
up such that thc small hole pockets around thc F point 
get filled up and the Fermi surface nesting disappears. 
Then the static "stripe" SDW order vanishes, and 
Dq will become unstable again. 

Rewriting 



ulu,^u}u,,ulu,,ul...ul,_^u, 



(16) 



with {is} ~ io, ii,..., iAT;^ a sequence of lattice sites con- 
necting i = iQ and j = i]\[.. and using Uj Ui^^^ ~ 
1 -f iBi^ ■ 77i,+i,i^ ~ e'''^^s■f^is+l.is (?7,^+i_,^ = r^^+i - r,J 
(neglecting the phase associated with the solid angle 
spanned by 11^^,11^^^^, z), one finds 



exp i dr T) 



exp 



exp 



dr ■ (DD) • dr' 



(17) 



where the spin correlation length ^ 
1/ [- J dujImV{0,iu)]^^^. Unless ImV{q,Lu) ^ 
at q — > as in the spin wave case, the spiral fiuctuations 
in (fT4|) will generally lead to a spin disordered state 
with a finite ^s, but its precise nature has to be self- 
consistcntly determined in view of the divergent I?(q, lu) 
in thc small q and uj regime. 
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The dynamic spiral fluctuations of the magnetization 
directions are kinetic energy driven by the electrons near 
the Fermi level. However, this will also make the elec- 
trons lose their long-wavelength coherence. Indeed, the 
single-particle propagator for the electrons to leading or- 
der of approximation can be written as 

G^(j,j;t)^([/t(t)[/,(0))G'^ (18) 

where is the propagator for the a-particles whose lead- 
ing term is coherent governed by the mean-field Hamilto- 
nian Hq in ([9]). Thus the dynamic spiral fluctuation will 
quickly damp the coherent motion of the quasiparticle 

beyond the spin correlation length via (^UJ{t)Uj{0)'^. 

By contrast, after averaging over the rii-field, the elec- 
tron singlet pair operator A^'" = J^a '^^la'^j-a '^^^ ^e 
expressed by 

cr 

« E'^«l°]-. (19) 

a 

at |i — j| ^ ^s- It clearly shows that in such a spin 
disordered state, the singlet pairs of electrons can still 
propagate coherently. The above contrast between the 
single electron and singlet pairs of electrons indicates an 
instability of the system towards superconductivity with 

Scr ^ ('^L'^]-ct) 7^ 0- The latter can be indeed reahzed 
by exchanging the dynamic spiral fluctuations, I?(q, a;), 
between the /3-electrons based on ([T2|) , which in turn sta- 
bilizes 2?(q, w) through a renormalized n^(q, cj) as the 
/3-electrons form Cooper pairs. It is noted that the sta- 
bihzed Rel?(q, w) < at small q and w will lead to an 
attractive (repulsive) interaction between /3-electrons if 
k + q is within (outside) the magnetic BZ according to 
([T2l) . It means that the /3-electrons will form the dom- 
inantly dj.2_,y2-wave pairing at four M points in Fig. 
2(b), similar to the case in the cuprates. Physically this 
kinetic-energy-driven pairing can be understood as that 
two electrons sitting at, say, NN sites, gain enhanced hop- 
ping energies by sharing the spiral twist of n,;'s between 
the two sites. 

Discussions. It is interesting to point out that Tc is 
upper bounded in the BCS theory because, whereas the 
softening phonon can enhance the attractive force, it also 
leads to the structural instability of the solids. Here 
the superconductivity is caused by strong dynamic spiral 
fluctuations which are in a "melting" SDW regime. Tc 
is expected to fall when such dynamic spiral fluctuations 
get reduced at higher electron doping concentration (the 
local magnetization Ma should be eventually destroyed 
when the Fermi level reaches beyond the highest /3-band 
shown in Fig. 1, presumably the one with the dominant 
d^2_y2 characters according to the LDA calculationfisj). 



Although the d-electrons in Fe-based compounds are 
believed quite itinerant, our proposal suggests that the 
underlying physics is still far from the conventional itin- 
erant magnetic metals in the following sense. Due to the 
presence of a large local SDW mean-fleld Ma, each band 
near the Fermi level are still split into the lower and up- 
per Hubbard bands, with the lower one filled by the elec- 
trons in both undoped and electron-doped cases, which 
are responsible for the origin of large Ma's. In contrast to 
the single-band Hubbard model relevant to the cuprate 
superconductors. [2^ however, here some of the upper 
Hubbard bands are already at the Fermi level even in the 
undoped case. A search for the depleted density of states 
below the Fermi level by, say, photoemission, [2^ may 
provide useful information concerning the correctness of 
the present model. The presence of a sizable local mo- 
ment, which is distinguished from the itinerant approach, 
may be also investigated via various magnetic measure- 
ment above Tc, including the magnetic susceptibility, [l^l 

Finally we caution that in the present model the effect 
of the unit cell doubling with two irons per cell due to the 
crystal field reason (As ions are displaced above and be- 
low the Fe plane alternatively) has not been considered, 
which may drastically affect our results if it becomes suf- 
ficiently strong because in that case the a- and /3-band 
splitting due to the local SDW may no longer be complete 
and Qf-bands can become partially filled. 

In conclusion, we have proposed a simple effective 
model to describe the low-energy physics in FcAs super- 
conductors. This model is based on an assumption that 
there exists a robust local SDW order, but its magneti- 
zation direction is intrinsically unstable against forming 
a true long range order due to the competition with the 
kinetic energy of the charge carriers near the Fermi en- 
ergy. In the undoped case, a weak SDW order of stripe 
type is formed due to the peculiar nesting structure of the 
Fermi surfaces. In the doped case, the dynamic melting 
of the SDW order will result in the d-wave superconduct- 
ing pairing of the doped electrons near the Fermi energy. 
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